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•  Latest  generation  of  high  power  battery  cells  are  often  comprised  of  multiple 
alternating  layers  of  anode,  cathode,  separator  and  electrolyte 

•  Macroscopically,  battery  cores  represent  an  orthotropic  material  subject  to  a 
time  variant  heat  source 

•  Safety  and  performance  considerations  place  a  premium  on  packaging  design 
and  installation  thermal  maintenance 


•  SAIC  has  developed  a  numerical  solver  tool  to  evaluate  thermal  performance 
of  packaging  alternatives  that  will: 

-  Run  rapidly  (avoiding  costly  finite  element  simulations) 

Evaluate  multiple  geometries  (prismatic,  cylindrical,  annular  cell  arrangements) 

-  Provide  flexibility  for  multiple  configurations  (air  or  liquid  cooling) 

-  Support  steady-state  and  transient  solutions 

-  Quantify  predictive  uncertainty 
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Allows  for  internal  cell  temperature  prediction  where  instrumentation  is 
difficult,  at  best 
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Cartesian  Steady-State  Solver 
for  Prismatic  Battery  Cells 


MODELING  mm  SIMULATION,  TESTING  m  VALIDATION 


m- 


•  Solver  predicts  temperature  response  within  the  battery  core  based  upon  user 
supplied  input  of  geometry,  cell  properties,  boundary  conditions  and  heating  rate 


User  can  select  mesh  refinement  and  degree  of  simulation  precision 
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If'  Cartesian  Solver 

for  Prismatic  Battery  Cells 


If 

.V 

MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


Solves  for  temperature  profiles  as  a  function  of 
coolant  temperature,  heating  rate,  effective 
boundary  conditions  and  cell  properties 
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The  solution  solves  for  temperature  response  with 
effective  convective  boundary  conditions  along  all 
six  faces  of  the  cell  core  utilizing  boundary 
conditions  of  the  third  kind 
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Similar  forms  exist  for  cylindrical  cells 
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|'A'  Cartesian  Steady-State  Solver 
for  Prismatic  Battery  Cells 


MSTV 

MODELING  m  SIMULATION,  TESTING  hud  VALIDATION 


Prismatic  Steady-State  Temperatures.vi 


Elle  Edit  Operate  Tools  Browse  window  tlelp 
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Prismatic  Cell  Steady-State  Core  Temperature  Response 


Temperature  Extrema 
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HDR  (W/C)  |  0.675 


—  Battery  Dimensions 
1 —  Battery  Core  Properties 
-  Effective  Heat  Transfer  Coefficients 
Problem  Specifications 

Number  of  Eigenvalues 
(per  coordinate  direction) 


Core  Temperature  Statistics 

Temperature  Extrema  and  Locations  — 1 

Energy  Balance  Solution  Check  — 1 

Dimensionless  Temperature  Profiles  Core  Face  Values  — 1 

T=f(x/L1,  y/L2)  Solution  Location  Monitor  — 1 

at  z/L3=(0.0,  0.1,  0.2,  0.3,  0.4,  0.5)  PlIOrTP 

- - GVScTS 


11  August  2011 


UNCLASSIFIED 


5 


^  Results  of  Several  Bas 
Verification  Simulation; 
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Uniform  heat  transfer  coefficients 

Note  increased  gradients  in  x-direction  due  to 
orthotropic  properties 


Adiabatic  conditions  along  x=0  face 


Adiabatic  conditions  along  y=0  face 


Adiabatic  conditions  along  x=0  &  y=0  face 

GVSETS 


UNCLASSIFIED 


6 


Solver  Example  -  Prismatic 
Cells  Between  Liquid  Coldplates 


MODELING  m  SIMULATION,  TESTING  min  VALIDATION 


•  A  generic  prismatic  cell  (15x10x1  cm)  with  sandwiched  between  cold  plates  - 
aluminum  interstitial  plates  act  as  cooling  fins  for  cell  lateral  surfaces 

•  Effective  heat  transfer  coefficients  estimated  through  analogous  thermal  circuits 

•  Similar  approach  used  for  y-direction  effective  heat  transfer  coefficients 


Monte  Carlo  Simulation  Tool 


MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


To  quantify  uncertainty,  a  Monte 
Carlo  simulation  tool  has  been 
developed 

User  selects  distribution  type 
and  parameters 

Currently  supports: 

-  Gaussian  distributions 

-  Log-normal  distributions 

-  Uniform  distributions 

User  tool  has  distribution 
samples  to  aid  user  identification 

User  selects  number  of  Monte 
Carlo  samples  to  simulate 


Monte  Carlo  Simulation 
Distribution  Parameters  Input  Screen 
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Monte  Carlo  Simulation  Results 


MODELING  mm  SIMULATION,  TESTING  m  VALIDATION 


Monte  Carlo  simulation  runs  multiple 
simulations  using  random  sampling  from 
defined  variable  distributions 

Simulations  run  rapidly  -  1000  samples 
runs  in  a  matter  of  minutes  on  a  laptop 

Allows  for  rapid  evaluation  of  uncertainty 

-  Material  properties 

-  Heating  rate 

-  Boundary  conditions 

Solver  plots  battery  cell  maximum, 
minimum  and  average  temperature 


Effective  HX  Coef  hi 


Effective  HX  Coefh2 


hi  (W/m2-K)  h2  (W/m2-K) 


Q  (W)  h3  (W/m2-K)  h4  (W/m2-K) 


Effective  HX  Coef  h5 


65  70  75  80  85  90 


Effective  HX  Coef  h6 


97.5  100  '  105  '  HO  '  115  '  120 


kperp  (W/m-K)  h5  (W/m2  K)  h6  (W/m2  K) 


34  36  38  40  42  2640  2660  2680  2700  272<J  1100  1150  1200  1250 

kpara  (W/m-K)  rho  (kg/m3)  Cp  (J/kg-K) 


Input  Property  Distributions 


Monte  Carlo  Simulation  Results 
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Monte  Carlo  Simulation  Results 


MODELING  fho  SIMULATION,  TESTING  m  VALIDATION 


•  The  solver  post-processes  the  Monte  Carlo  simulation  results  to  give: 

Distributions  for  the  minimum,  average  and  maximum  core  temperatures 

-  Heat  Dissipation  Rate  (HDR)  distribution 

-  Parameters  of  those  distributions  (mean  and  standard  deviation) 


Tmax  mean 

69.77 

Tmax  std  dev 

1.83 

Tave  mean 

68.03 

Tave  std  dev 

1.80 

Tmin  mean 

64.39 

Tmin  std  dev 

1.72 

HDR  mean 

0.67 

HDR  std  dev 

0.04 

Distribution 

Parameters 


Heat  Dissipation  Rate  (HDR)  is  a  measure  of  the  package  cooling  effectiveness 
and  is  defined  as:  Q 


HDR  - 


T  —T 

Max  Coolant 
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•  Heat  Dissipation  Rate  (HDR)  can  be  used  to  quantify  the  expected  temperature 
difference  from  the  cell  to  the  coolant  as  a  function  of  heating  rate 

•  HDR  leads  to  predictive  of  maximum 
and  minimum  core  temperatures  as  a 
function  of  heating  rate 

with  uncertainty  bounds  (±3 a)  TemDeratur6 


Temperature 
Difference 
to  Coolant 


HDR  = 


T  —T 

Max  Coolant 


•  This  illustrates  best-  and  worst-case  Ceii  Heating  Rate  (v 

cooling  scenarios  for  this  particular  packaging  design 

•  For  example,  at  a  30°C  coolant  temperature  and  a  15W  cell  heat  load 

-  ATmax  (Cell  max  temperature  -  Coolant  Temperature)  =  27  C 

Worst-case  maximum  core  temperature  is  predicted  to  be  48  <  Tmax  <  57°C. 
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Cartesian  Transient  Solver  for 
Prismatic  Battery  Cells 


>  Prismatic  Transient  Tenipcraturcs.vi 
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Prismatic  Cell  Transient  Core  Temperature  Response 


Z-plane  Temperature  Contours  (0.0, 0.26, 0.60) 
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•  Allows  for  transient  simulation  of  time-variant  heat  loads 

•  Heat  loads  can  be  user  defined  or  uploaded  from  tab-delimited  files 
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Transient  Results 
Step  Discharge  Case 


MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


iw 


•  Results  of  a  step  change  in  heating  from  initial  conditions 


Battery  Core  Temperature  Response 


•  Identifies  battery  packaging  concept  time  constant,  worst  case  loading 
expectations  and  recovery  time 
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gp  Transient  Results 
f/’  Load  Leveling  Case 


MODELING  fho  SIMULATIDN,  TESTING  m  VALIDATION 


•  Load-leveling  presents  a  challenging  thermal  demand 


•  This  case  represents  the  temperature  response  to  repetitive  cycles  of  60  sec 
discharge  (100A)  and  120  sec  charge  (50A) 


Battery  Core  Temperature  Response 
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Transient  Results 
Driving  Scenario  Case 


MODELING  fho  SIMULATION,  TESTING  m  VALIDATION 
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•  Realistic  driving  scenarios  imposed  on  a  battery  used  for  mobility  assist 


•  Profiles  typically  show  extensive  non-uniformity 


100  200  300  400  500  600  700  800  900  1000  1100  1200  1300  1400 

Time  (sec) 

Battery  Core  Temperature  Response 
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Wrap  Up 


MODELING  HMD  SIMULATION,  TESTING  hud  VALIDATION 


•  Thermal  solver  predicts  internal  battery  core  temperature  response 

•  Steady-state  and  transient  applications 

•  Different  geometries  supported  -  prismatic,  cylindrical,  annular 

•  Monte  Carlo  simulator  included  to  quantify  uncertainty 

•  Several  future  developments  are  envisioned: 

-  Link  to  battery  performance  tool  (electrical  model) 

-  User  tools  to  support  boundary  condition  estimation 

-  Inclusion  of  transient  coolant  temperatures 

-  Temperature  dependent  property  effects 

-  External  package  thermal  inertia  effects 

•  Questions  /  Comments  /  Feedback 
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